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We study the effect of a one dimensional optical superlattice on the superfluid fraction, number 
squeezing, dynamic structure factor and the quasi-momentum distribution of the Mott-insulator. 
We show that due to the secondary lattice,there is a decrease in the superfluid fraction and the 
number fluctuation. The dynamic structure factor which can be measured by Bragg spectroscopy 
is also suppressed due to the addition of the secondary lattice. The visibility of the interference 
pattern (the quasi-momentum distribution) of the Mott-insulator is found to decrease due to the 
presence of the secondary lattice. Our results have important implications in atom interferometry 
and quantum computation in optical lattices. 
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I. INTRODUCTION 



When a gas of ultracold atoms is loaded into an optical lattice, its properties are modified strongly [ill. Ultracold 
bosons trapped in such periodic potentials have been widely used recently as a model system for the study of some 
fundamental concepts of quantum physics like Josephson effects[l|, squeezed states, Q landau-Zener tunneling and 
Bloch oscillations 3 and superfluid-Mott insulator transition Using superposition of optical lattices with different 
periods it is now possible to generate more sophisticated periodic potentials characterized by a richer spatial 
modulation, the so-called optical superlattices. An important and exciting application of optical superlattice is 
quantum computation 0]. The physics of one-dimensional optical superlattices has been a subject of recent research, 
including fractional filling Mott insulator (MI) domains [1], dark Qand gap fiol solitons, the Mott-Peierls transition 
11 1 , non-mean field effects 1131 , phase-diagram in two colour superlattices [IjI , Bloch-Zener and dipole oscillations 
14 1, collective oscillations jl5l|and Bloch and Bogoluibov spectrum [l^. A key observable in these systems is the 
interference pattern observed after releasing the gas from the lattice and letting it expand for a certain time of flight. 
Monitoring the evolution of this interference pattern reveals e g., the superfluid fraction, number squeezed states 
quasi-momentum distribution, observation of collapse and revivals of coherence due to atomic coherence [l^ 
and superfluid to Mott insulator transition [l^ . Further coherence properties of Bose-Einstein condensates offer 
the potential for improved interferometric phase contrast. The MI state plays a central role for various quantum 
information processing schemes j2^]. Because of the experimental importance of BEC in optical lattices, it is crucial 
to understand the influence of the secondary lattice which is emerging as a new manipulating tool on the coherence 
properties of a BEC. In the present paper, we study in what way the superfluid fraction, number fluctuation, the 
dynamic structure factor and the quasi-momentum distribution (and hence the visibility of the interference pattern) 
of the MI is influenced by the addition of the secondary lattice. 



II. THE BOGOLUIBOV APPROXIMATION TO THE BOSE-HUBBARD HAMILTONIAN 

The light shifted potential of the superlattice is described as 



V{z) = V, cos2 + V2 cos' (^^ + (1) 

Here di and ^2 are respectively, the primary and secondary lattice constants. Vi and V2 are the respective ampli- 
tudes. The secondary lattice acts as a perturbation and hence we will take V2 << Vi. is the phase of the secondary 
lattice. When cf) — 0, each site of the lattice is perfectly equivalent due to the symmetries of the system so that the 
population and on site energies are same at each site. An asymmetry is introduced when (j) ^ Q and hence the onsite 
energies are not the same at each site. The harmonic trapping potential is given by Vho{r, z) — ^ (uj'^r'' + co'^z') and 
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the optical superlattice potential is given as Vop = En (si cos^(^) + S2 cos^(H). In our case we take ^2 = 2di = 2d 
which gives rise to a periodic double well potential. Also si and S2 are the dimensionless amplitudes of the primary 
and secondary superlattice potentials with si > 82- Er = is the recoil energy (wfl = is the correspond- 

ing recoil frequency) of the primary lattice. [/ — is the strength of the two body interaction and a is the two 

body scattering length. We take ojr > oJz so that an elongated cigar shaped BEC is formed. The harmonic oscillator 

frequency corresponding to small motion about the minima of the optical superlattice is ojs ~ md'^ ■ ^^'^ BEC 
is initially loaded into the primary lattice and the secondary lattice is switched on slowly so that the BEC stays in 
the vibrational ground state. The frequency of each minima of the primary lattice is not perturbed significantly by 
the addition of the secondary lattice. lUs » uJz so that the optical lattice dominates the harmonic potential along 
the z-direction and hence the harmonic potential is neglected. Moreover we also take a sufficiently large harmonic 
confinement in the xy plane which effectively reduces the problem to one-dimension. The strong laser intensity will 
give rise to an array of several quasi-two dimensional pancake shaped condensates. Because of the quantum tunneling, 
the overlap between the wavefunctions between two consecutive layers can be sufficient to ensure full coherence. Fol- 
lowing our earlier work fl^ the effective one-dimensional Bose-Hubbard Hamiltonian for I lattice sites and 7^ is 
written as 



H 



U' 



E 
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Here Jj is the site dependent strength of the Joscphson coupling and is different when going from j — 1 to 
j and j to j + l.The two Josephson coupling parameters are conveniently written as Jq ± Aq/2, where Jq — 

^ y/si ~ si exp ^— ^ and Ao = S2Eftexp ^— ^^'^ ^ . The strength of the effective on-site interac- 
tion energy is f/e// = U J dz \w{z)\'^. Here C^e// = C4///V2d, V2d is the two dimensional area of radial confinement (i.e 
area of each pan cake shaped BEC). €j is the onsite energies and takes two distinct values (ei and €2) corresponding 

to odd and even sites. In the mean field approximation, the operators dj and aj are classical c numbers, dj = 
Stationary states with a fixed total number of particles N are obtained by requiring that the variation of iJ — 
with respect to (p* vanish. Here ji is the chemical potential. This yields the eigenvalues equation 

+ Ueff 1^ (pj - Jjpj+i - - fJ-Pj = 0. (3) 

We write (pj as 

= gje'''""- (4) 

The eigenvalues are found as 

2Ueffno - ^[Ae]' + 4el 
A^- ^ . (5) 

Where = \J 4 Jq cos^ 2kd + Aq sin^ 2kd and Ae = £1 — 62. The eigenvalue fi corresponds to the chemical potential 
for k = 0. The Bogoliubov spectrum of elementary excitation describes the energy of small perturbations with quasi- 
momentum q on top of a macroscopically populated state with quasi-momentum k. In the Bogoliubov approximation, 
we write the annihilation operator in terms of c-number part and a fluctuation operator as 

dj = (0 + Sj) exp f ) (6) 



The resulting Bogoliubov equations for the fluctuation operator Sjin the optical superlattice take the following form 

ihSj = {2Ueffnn - ^) Sj - JjSj+i - J,-_i(5j_i + UeffriQ^ (7) 
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Figure 1: The acoustical branch of the Bogohubov spectrum as a function of qd. naUef f / Er = l,Jo/Eu = 1,Ao/Er = 0.1. 
The bold curve is for (ei — e2)/ER = 1, while the dashed curve is for ei = e2.The change in the Bogoliubov spectrum due to 
asymmetry is appreciable only near the band edge. 

no is the 2d average density of atoms per site of the lattice. The above equation is solved by constructing quasi- 
particles for the lattice, which diagonalizc the Hamiltonian i.e 



(8) 



The quasi-particles obey the usual Bose-commutation relations 



The excitation amplitudes obey the periodic boundary conditions 



(9) 



q q q q 

Finally the phonon excitation frequencies are found to be 

1 



Where 



Ae 1 



/3i = Ueffuo - — + ^^{Aef + 16J2. 



(10) 



(11) 



(12) 



h = Ueffno + ^ + i^(Ae)2 + 16j2. 



(13) 



= Y 4Jo cos2 2qd + Ag sin"' 2qd 



(14) 



where fUVq,- is the acoustical branch. There is another branch called the gapped branch (analogue of the optical 
branch) whose energy is given by hu}q,+ . In Figure [1] we find that for Ae not large, the asymmetry due to 7^ 
does not appreciably change the Bogoluibov spectrum.lt is only near the edge of the Brilliouin zone that a small 
change is visible. 



4 



Due to the above fact in the following we will only discuss the simple case when (/> = i.e experimentally <j) does 
not deviate much from zero so that ei =62- This case also allows us to tackle the problem analytically. The resulting 
equations for amplitudes and frequencies are solved to yield the Bogoliubov amplitudes as 



^ I 2Jocos2qd + iAosm2qd \ ^^^^p ^^^^ 



y/ 4 Jq cos 2qd 


+ A§ sin 2qd 


2Jq cos 2qd + 


i Aq sin 2qd 



^ ~ 1^ V4J^os2<zd + Ai^sin2gd ) ' ^ ' ^ ^ 



go go 



(19) 



where = lq-(2nQU f + Cg,-) and e^^- = 2 Jq — y 4 Jq cos^ 2qd + Aq sin 2qd. 

III. SUPERFLUID FRACTION AND NUMBER FLUCTUATIONS 

An interacting many body system is said to be superfluid, if a condensate exists. This happens when the one-body 
density matrix has exactly one macroscopic eigenvalue, which defines the number of particles in the condensate. The 
corresponding eigenvector describes the condensate wavefunction, "iii^if) = e*"^''"^|V'o(^)P- The superfluid velocity is 
given as 

Vs(r) = -^V0(r) (20) 
m 

Here m* is the effective mass of a single atom in the optical superlattice. We now write down the expression for 
the superfluid fraction based on the rigidity of the system under a twist of the condensate phase [2lj | . Suppose we 
impose a linear phase twist (/)(r) = with a total twist angle Q over a length L of the system (with ground state 
energy Eq) in the z direction. The resulting ground state energy, Eq will depend on the phase twist. Thus, 

Ee-Eo = ^m*7V/,«2 (21) 

where N is the total number of particles, fs is the superfluid fraction and m* = 2d'^{4ji-A'^) ■ Substituting equation 
(21) into (22)gives 

_ AJojEe - Eq) 

7V(4J2- A2)(A0)2 ^^^^ 

Here A9 is the phase variation over 2d. We now need to calculate the energy change {Eq — Eq) using second order 
perturbation theory, under the assumption that the phase change, A6 is small. This yields 

{Eg - £^0) = AE^^^ + AS(2) (23) 



Where AE^^^ is the first order contribution to the energy change 



Here \tpo) is the ground state of the Bose-Hubbard Hamiltonian. The hopping operator f is given by 



1 



The second order contribution is written as 



where the current operator J is 



1 



The total superfluid fraction has two contributions. 



where 



(24) 



(25) 



(26) 



(27) 



(28) 



2Jo 



N{Ui-Al) 



V-olTIVo 



(29) 



72 _ 2^ 



N{4J§ - Ag) - Eo 



Using the expressions for the various Bogoliubov amplitudes and frequencies, we can now evaluate /i^^ and /. 



(30) 



p(2) 



2Jn 



(31) 



In the Bogoliubov approximation this takes the form 



^^'^ = jV(4Jg - Ag) ^ + ^^+^^^- + ^^^^+^1^° 

The fluctuation operators appearing in equation (32) are now written in terms of the quasi-particle operators. 



(32) 



(33) 
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Figure 2: The superfluid fraction as a function of S2/S1 with noUef f / Jo = 1 , / = 10 and n = 10. As the strength of the 
secondary lattice increases with a fixed strength of the primary lattice, there is a quantum depletion of the condensate which 
is seen as a decrease in the superfiuid fraction. 



Finally, we find in the zero temperature limit 



mUf^) + E ^0 {ulu.e^"^'^ + u,ule--^^) | (34) 

Here, the summation runs over all quasi- momenta q = jji with j = 1, 2, ...(/ — 1). The normalization condition is 
obtained by putting /i^'' = 1 when d 0. 



Y: JM) + JoY: Jo2Reiu,u;) = ^^'f^ (35) 

Using the Bogoliubov amplitudes derived in the previous section, one can show that fs — 0. Consequently, we 
find that the total superfluid fraction has contribution from just fs^K A plot (Figure [2]) of the superfiuid fraction as a 
function of S2/S1 reveals a decrease in the superfiuid fraction as the strength of the secondary lattice increases. This 
is to be expected since in the presence of the secondary lattice, it has been shown that there exists a fractional filling 
Mott insulating state in the phase diagram [8j].This itself is an indication of a reduced superfiuid fraction. This result is 
in accordance with earlier work of Rousseau et al. [l^l where they have considered the effect of a secondary lattice on 
an one dimensional hard core bosons(strongly correlated regime). As the strength of the secondary lattice increases, 
we approach the Mott-insulator transition. Since the phase twist is equivalent to the imposition of an acceleration 
on the lattice for a finite time, the condensate now in the superlattice seems to resist this acceleration or simply tries 
to resist the phase twist and thus there is a reduction in the superfiuid flow. A direct consequence of the decrease 
of the superfluid fraction is a decrease in the number fluctuation, which we show below. Increasing the lattice depth 
reduces the tunneling rate between adjacent wells. This can be viewed as a reduction of the number fiuctuations at 
each lattice site. As the probability of the atoms to hop between wells decreases, the number variance tT„goes down. 
Quantum mechanically, this implies that the phase variance describing the spread in relative phases between the 
lattice wells, has to increase. This effect can be seen directly by looking at the interference pattern of a BEC released 
from an optical trap. We can find an expression for the fluctuations in the relative number in each well as [2l| 
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(36) 



and 
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Figure 3: The number fluctuation as a function of S2/S1 with nof/e///Jo = 1 , / = 10 and 71 = 10. As the strength of the 
secondary lattice increases, there is a loss of superfluidity. The interplay of the interaction and tunneling terms renders number 
fluctuations energetically unfavorable. The number fluctuations decrease with increasing potential of the secondary lattice. 
There is a corresponding increase in the phase fluctuations. 



{Ug 



(37) 



/ is the total number of sites and n is the mean number of atoms on each site of the lattice. A plot (Figure [3])of 
the number fluctuations versus S2/S1 reveals as expected a decrease with increasing strength of the secondary lattice 
indicating a loss of phase coherence. The number variance may be measured experimentally by studying the collapse 
tc and revival trev times of the relative phase between sites [23| . The relation is given by cr„ = |^ . This reduction 
in the number fluctuation is also called as the atom number squeezing. This increased squeezing as a result of the 
secondary lattice has an important application in in improved atom interferometry since with increased squeezing the 
coherence time also increases [IJ]- These atom number squeezed states have reduced sensitivity to mean- field decay 
mechanisms. The secondary lattice then serves to coherently maintain a balance between coherence as well as the 
decoherence effects due to mean-field interaction. 



IV. DYNAMIC STRUCTURE FACTOR 



The capability of the system to respond to an excitation probe transferring momentum p and energy huj is described 
by the dynamic structure factor. In the presence of a periodic potential the dynamic structure factor takes the form 



S{p,Lo) = Zo.{p)S[lu - uj^ip)] (38) 

a 

where Za{p) are the excitation strengths relative to the a*'* mode. a is the band label. For each value of the 
quasi-momentum q, there are infinite set of excitation energies hua{q). It is often convenient to consider values of 
q outside the first Brillouin zone and to treat the energy spectrum and Bogoliubov excitation amplitudes ^ and 

a periodic with period 2qB- Here 9b = is the Bragg momentum denoting the boundary of the first Brillouin 
zone, p is assumed to be along the optical lattice (z axis), is not restricted to the first Brillouin zone since it is the 
momentum transferred by the external probe. The quantities q, p and qs are related as q = p + 2lqB, I is an integer. 
In the first Brillouin zone I = 0. The excitation energies huJa(j>) are periodic as a function of p but this is not true for 
the excitation strengths Z^. The excitation strengths can be evaluated using the standard prescription [25j 



Zaip) = 



(39) 
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Since l^j „P=|w'_|_2_(jP and |v'„P=|?^J_|_i j^P, we will drop all j dependence from the Bogoliubov amplitudes. The 
excitation frequencies for different a has already been derived in our earlier work .[l^ We are interested in the low 
energy region where Zi{p) is the dominating term arising from the first band. The dispersion law for the lowest band 



huJiip) = y^ep(2not/e// + ?p) (40) 



ep = 2 Jo - ^4J2 cos2 (^^) + Al sin^ (^^) (41) 

The behaviour of Zi{p) can be studies analytically in the tight binding limit. In this limit one can approximate the 
Bogoliubov amplitudes in the lowest mode as. 

u„(z)=^e'-'W/(^-2j-rf) (42) 
j 

and analogously for Va{z), where f{z) is a function localized near the bottom of the optical potential V at z — 0, 
and J labels the potential wells. Within this approximation the function / also characterizes the ground state 
order parameter which reads 4>{z) — f{z — 2jd). We can approximate the function f(z) with the gaussian 

f{z) = exp [— z^/2cr^] / (n^^'^y/a). The width a is found by minimizing the ground state energy 



and behaves like a ^ -, — ■ — '',,,1 /. . After some trivial algebra we find 



dz (43) 



The expression for Z-y (p) shows both the oscillatory behaviour through f^J^^^p-^ and decaying behaviour at large p 

/ 2 2 2 \ 

through exp ( — gj^^^f j ■ Figure [4] shows the excitation strength Zi{p) for two values of |^ = 0.1 (solid line) and 

^ = 0.4 (dashed line). On increasing the strength of the secondary lattice, Zi{p) is quenched. This behaviour can 

be understood by looking at the low p limit of S{p) = J S{p,uj)du} = J^^ [/ . • '-''^ increasing S2, m* increases 

and hence S{p) decreases. The presence of the secondary lattice results in the suppression of Zi{p). The system now 
becomes more heavy and is not able to respond to an external excitation probe. The momentum transferred is now 

comparatively less. Note that in the absence of interactions, the oscillatory behaviour disappears and the strength 

/ 2 2 2 \ 

reduces to Zi(p) = exp ( — ^^2^^ j ■ This shows that the effect of the secondary lattice on the quenching is present only 

in the presence of interactions. The zeroes of Zi (p) a,t p — 2lqB reflects the phonon behaviour of the excitation spectrum 
which also vanishes at the same values. The quantity Zi(p) can be measured in Bragg spectroscopy experiments by 
applying an additional moving optical potential in the form of Vsit) — Vbcos(^) — cot. The momentum and the 
energy transferred by the Bragg pulse must be tuned to the values of p and hoj corresponding to the first Bogoliubov 
band. 

V. QUASIMOMENTUM DISTRIBUTION OF THE MOTT INSULATOR IN AN OPTICAL 

SUPERLATTICE: VISIBILITY OF FRINGES 

For a Bose-Einstein condensate released from an optical lattice, the density distribution after expansion shows a 
sharp interference pattern. In a perfect Mott-insulator, where atomic interactions pin the density to precisely an 
integer number of atoms per site, phase coherence is completely lost and no interference pattern is expected. The 
transition between these two limiting cases happens continuously as the lattice depth is increased. In this section, 
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Figure 4: The excitation strength Zi{p) for two values of |^ = 0.1 (solid line) and ^ — 0.4 (dashed line).l7e///Jo ~ 0.2. 
The figure shows both the oscillatory behaviour through i^^^^p-^ and decaying behaviour at large p through exp ^— '^^2^-? j.On 

increasing the strength of the secondary lattice, Zi{p) is found to be quenched. The first maxima is found near the edge of the 
first Brillouin zone. 



we will look into the influence of increasing the strength of the secondary lattice on the phase coherence of the 
insulating phase. We consider an integer number n of atoms per site and Jq ± ^ << Ueff- In this situation the 
gas is in the Mott-insulator phase. The Mott insulating phase has the property that the fluctuations in the average 
number of particles per site goes to zero at zero temperature. These fluctuations can be described as quasihole and 
quasiparticle excitations. To calculate the quasimomentum distribution S{k) for a flnite tunneling, path integral 
techniques can be applied to obtain the single-particle Green function, G{k,uj). The quasi-momentum distribution is 
an useful quantity to describe the interference pattern observed after release of the cold cloud from the optical lattice. 
From the absorption image of such an interference pattern, the phase coherence of the atomic sample can be directly 
probed. To extract quantitative information from time-of-flight absorption images, one can use the usual definition of 
the visibility of interference fringes [1^ , 



= s — 

The quasimomentum distribution iS'(fc) contains information about the many-body system which is periodic with 
the periodicity of the reciprocal lattice corresponding to the secondary lattice. Thus to predict the interference pattern 
in the superlattice, our goal is to calculate S{k) as function of Jq and Aq. We calculate the quasiparticle and quasihole 
dispersions using the functional integral formalism of Van Oosten et. al. (27j . The grand-canonical partition function 
in terms of the complex functions a* (t) and Oj (r) is written as 



Z = Tre-^^"^ J Da* Da exp {-S [a*, a]/ h} (46) 
where the action S[a* , a] is given by 



.hp 

S[a* , a] — dr 
Jo 



(47) 



Jjj' is the hopping element, (3 = l/fcsT, fc^ is the Boltzmann constant and T is the temperature. A Hubbard- 
Stratonovich transformation decouples the hopping term. 



S [a* , a,ip* ,ip] — S [a* , a] + 



hp 



dr 



3,3' 



(48) 
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Here ip* and V are the order parameter fields. Integrating over the original fields Uj and aj, we find 

exp{-S''ff[il)*,i}]/h) = exp i-\J^ ^ dTY^Jjj,i>*i>j, \ J Da* Daexp(^-S^°^[a*,a]/h^ 



exp 



1 f'^'^ i 

J,]^ I " 5 ^ 



(49) 



Here S^^^[a*,a] is the action for Jjj' = 0. We can now calculate S'^^f perturbatively by Taylor expanding the 
exponent in the integrand of equation (49) and find the quadratic part of the effective action using {aj(ij')g(o) — 

S(^\rM = j drl Y,rj{r)^rir) --j dr' Jjj'JwrAr) i %(r)<(r'))s(o, Vi'(r') (50) 



We first evaluate the part linear in J^y for nearest neighbours. We have 

^ Vr(r)^,'(r) = [,h + ^) V'.Vi+i + U - ^ ) rii>3-r 



(51) 



We now introduce V'j = [ufc + i(— l)-'wfc]ea;p(ij2fc(i). As the condensate moves from one well to the next, it acquires 
an additional phase, which depends on the height of the barrier. As the height alternates and hence the tunneling 
parameter, the phase also alternates. This picture is conveniently represented by the j dependent amplitude. This 
implies 

^■(/)j*(T)V'j'(r) = 2 Jo - |ufc|^] cos(2A;d) - VXJq [uy^v]. + u^Ufc] cos(2fcc!) +iAo [|wfe|^ - \vk^\ sin(2fcc!) 
i.i' 

+Ao {ukvl + UfcUfe] sin(2fcrf) (52) 
For the imaginary part to vanish we have for the one-dimensional optical lattice 

,* , Aosin(2A;d) 

u^Vk = UkVk = i^ki^k ^ (53) 



2Ao cos{2kd) 



ek 



(54) 



efe = J4J§ cos2(2fcd) + A§ sin^(2A:d) 



Finally we have. 



(55) 



Next we calculate the part that is quadratic in Jj^ji . We can treat this part by looking at double jumps. 

JjyJii'ip*,{T) (aj(r)a*(r'))s(o) V'i'(^') = ("j Ws*(^')>s(o) Y •^j/'^Ji'V'l' WV'i'(T') 



j'i' 



(56) 



j'j' 



(57) 
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The first term in the summation is a jump forward, followed by a jump backward. The second is two jumps in the 
same direction. The above quadratic term then reduces to 

J,j'Ju'^*,iT) K(r)a*(T'))5(o, A'ir') = (a,(r)a*(r'))5<o, ^ eI^:(r)Vfc(r') (58) 

j'ii' k 

The Green's function is then easily calculated by following the steps indicated in ref.[27| 

G{k,Lu) _ Zk I - Zk 



hu + fi- E^' kuj + ^- El 



(59) 



The quasiparticle energies E^ are derived as 



= -| + Ueff - ^) ± ^ W 4 - 4e.t/e// (" + ^) + Ulff (60) 



The particle weight Zk is 



Zk - I , (61) 

The quasimomentum distribution can be directly calculated from the Green function G{k^uj) using the relation 

— G{k,Lo)exp{-iuj5t) (62) 
27r 

This yields 



S{k) is simply the quasi-momentum distribution which tells us about the many-body system. The visibility of 
the interference pattern of a cloud of BEG released from an optical superlattice as a function of the strength of 
the secondary lattice is shown in figure 4. As the strength of the secondary lattices increases, the visibility worsens 
since the system gradually goes deeper into the Mott insulator regime and a corresponding gradual loss of long range 
coherence. A finite visibility even for a Mott-insulator is due to short range coherence since the system consists of 
a small admixture of particle-hole pairs on top of a perfect Mott-insulator. A loss of visibility in the superlattice 
naturally means that there is loss of particle-hole pairs. 



VI. CONCLUSIONS 



We have studied the effect of a one dimensional optical superlattice on the superfluid fraction, number squeezing, 
dynamic structure factor and the quasi-momentum distribution of the Mott-insulator. We have shown that the 
secondary lattice suppresses the superfiuidity due to quantum depletion of the condensate and hence generates atom- 
number squeezed state which offers a possibility to create states with reduced sensitivity to mean field decay mechanism 
useful for improved atom-interferometry. A coherent control over the phase coherence in the superfluid as well as the 
Mott-insulating state can be achieved which has important applications in quantum computing. 
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Figure 5: The visibility of the interference pattern produced by an ultracold cloud released from an optical superlattice as a 
function of S2/S1 with Ueff/Jo — 40 and no = 3. As the strength of the secondary lattices increases, the visibility worsens since 
the system gradually goes deeper into the Mott insulator regime and a corresponding gradual loss of long range coherence. 
A finite visibility even for a Mott-insulator is due to short range coherence since the system consists of a small admixture of 
particle-hole pairs on top of a perfect Mott-insulator. A loss of visibility in the superlattice naturally means that there is loss 
of particle-hole pairs. 
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